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Abstract 

We study the static dielectric function e(q) of 2D electron system with 
spin-orbit coupling in the frame of the random phase approximation. We 
demonstrate that, in addition to the well-known 2/cp-Kohn anomaly, spin- 
orbit coupling gives rise to the novel anomaly in the dielectric function at 
small q = qo <C k>F, where qo is the distance between two Fermi surfaces. As a 
result of this anomaly a large-distance behavior of the potential from a point 
charge exhibits (in addition to the conventional Friedel oscillations) novel os- 
cillations with a period 2ir /qo . The remarkable feature of these oscillations is 
that they are not smeared out by the temperature. We show that the small-g 
anomaly also modifies the indirect exchange interaction of localized magnetic 
moments (RKKY interaction). In the presense of spin-orbit coupling this 

interaction acquires a high-temperature component. 
PACS number(s): 71.45. Gm,73.20-r, 73.20.Dx 
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I. INTRODUCTION 



The dielectric function of a two-dimentional electron system was first evaluated by F. 
Stern [1] within the random phase approximation (RPA). In the static case, Stern obtained 
the following expression for e(q) 

e(q)=e (l-v(q)F(q)), (1) 

where q is the wave vector and v(q) = 2-ne 2 je^q is the Fourier component of the Coulomb 
interaction, So is the dielectric constant of the medium. The response function F(q) is given 
by 

F{q) = -v (q< 2k F ) (2) 



F(q) = -?j(q-y/q*-4k* F ) (q>2k F ), (3) 

where v = m* /nh 2 is the 2D density of states; m* is the effective mass. It is well-kown that 
the singularity in F(q) at q = 2k F leads to the Friedel oscillations of the screening potential. 
For a point charge at a distance d from a 2D plane the oscillatory part of the potential has 
the form 

U-) = -(- ) sj £r4^* d , (4) 

\e a B J (2k F p) 2 

where as = h 2 eo/m*e 2 is the effective Bohr radius. Eq. (§) applies at large distances, 
k F p ^> 1. The result (|j) is derived for zero temperature T = 0. Obviously, the oscillations 
are smeared out as T increases. This effect amounts to an additional factor A(u) in Eq. 
defined as 

a 

A(u) = —— v (5) 
smh(w) 

where u = nk F pT/ E F . It follows from Eq. (5) that the power-law decay of V osc (p) persists 
at finite temperature from p ~ kp 1 to p ~ k F (E F /T). At larger p a crossover to the 
exponential decay cx exp(—2nk F pT/E F ) takes place. 



The expression for V osc (p) can be straightforwardly generalized to the case of two or more 
2D subbands. This is because the size- quantization wave functions in the direction perpen- 
dicular to the plane are orthogonal, so that the inter-subband transitions are forbidden. In 
this case V osc (p) is described by the sum of the Friedel oscillations from each subband. 

The same argument would apply in the presence of the Zeeman splitting if one neglects 
the orbital action of the magnetic field. In this case, the transitions between the branches 
of the spectrum with different spin projections are forbidden. 

In the present paper we demonstrate that the situation is quite different in the presence 
of the spin-orbit (SO) coupling. With SO coupling, the electronic states are classified not 
by the spin orientation but by the chirality. Our main results are as follows: (i) SO coupling 
gives rise to a new type of oscillations in V(p) which are caused by virtual transitions between 
the states with different chirality; (ii) the distinguishing feature of these oscillations is that 
they survive the increasing of temperature (not smeared out with the increasing T). 

To account for SO coupling, we choose the simplest form of the SO interaction for 2D 
electron system [2] 



where k is the momentum, n is the unit vector normal to the 2D plane, a = (a"i, cr 2 , 03) are 
the Pauli matrices, and a is the SO coupling constant. With spin-orbit term Eq. (6) the 
splitted Fermi surfaces for a given 2D subband represents two circles with radii (see Fig. 1) 



H S o = ctk ■ (cr x n), 



(6) 




(7) 




(8) 



where fop is related to the Fermi energy as 




(9) 



and q = k F ' — k F ' is expressed through the SO coupling constant as follows 
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2m* 



a 



(10) 



go = 



ti 
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It is important that the difference k F — k F does not depend on the position of the Fermi 
level. 

Our main observation is that, in contrast to the case of different 2D subbands, the virtual 
transitions between SO-splitted states within the same subband are allowed and, thus, do 
contribute to the response function. Below we will demonstrate that this contribution is 
singular at q = q . This singularity results in the following oscillatory behavior of the 
screening potential 



Comparing ([TT]) to @, one can see that the SO-induced oscillations do not contain the factor 
exp(—2kpd). Strictly speaking, this factor should be replaced by exp(-qod) in Eq. (fLl]) , but 
we assume throughout this paper that SO coupling is weak, i.e. q -C k F . The magnitudes 



elevated temperatures, Vso will always dominate since the conventional Friedel oscillations 
get damped. 

It is interesting to note that as far as the conventional Friedel oscillations are concerned, 
there is also a difference between the case of two size-quantization and two SO-splitted 
subbands. In the former case the large-p behavior of the screened potential is V(p) oc 
(sm(2kp p) + sin(2kp^ p)) / p 2 , whereas in the latter case, as we will demonstrate below, these 
oscillations have the form V(p) oc (sin^? + kf)p)/p 2 . 

The paper is organized as follows: In Sec. II we evaluate the dielectric function e(q) for 
different domains of wave vector q. The spatial oscillations of the screened potential are 
studied in Sec. III. Concluding remarks are given in the Sec. IV. 




(11) 



of the oscillations Eq. (|TT| ) and Eq. (^) are of the same order for d = 0. However, at 



II. DIELECTRIC FUNCTION WITH SO COUPLING 



4 



A. Spectrum and Eigenfunctions 



With SO term Eq. (IB), the Hamiltonian of a free electron has the form 



H = h ak ■ {<t x n). 

2m v ' 



The energy spectrum and the wave functions are given by 



2 V.2 



2m 



jj,ak 



(12) 



(13) 



V2 



i\xe 



(14) 



Here = ±1 is the chirality, A is the normalization area, and 0k is the azimuthal angle 
of the vector k. 



B. General Formula for the Dielectric Function 

We start with the RPA formula for e(q) applied to the case when free electron wave 
functions are the eigenfuctions of the Hamiltonian Eq. (|i~2"D 

(15) 



e M = 1 - v(a ) V Ivt v ,2 "(^(fc))-n(^(|k + q|)) 
, " k^ " fy(|k + q|) ' 

where n(E^(k)) is the Fermi-Dirac distribution function, and 



IXmXicW I 2 = + /V cos (0k+ q - 0k)), (16) 

are the overlap integrals. It is important to note that the virtual transitions between different 
chiralities are allowed. In the case of two size-quantization subbands, the overlap integral 
would reduce to the Kronecker symbol, so that e(q) would represent a sum of contributions 
from each subband. In contrast, as it follows from Eq. fll6|), the overlap integral vanishes 
for jj, ^ jj! only if k || q. 

In evaluating the dielectric function, it is convenient first to perform the angular inte- 
gration by introducing the auxilary variable in Eq. (|i~5|) 
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i-^y [°° dkk Hdpp 

^ lo h Jo 



Att 2 



^ / ,k-(k + q) 
dM l + fjifi — ^ — I x 
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kp 



n(^(fe))-7i(^(p)) \ , 



where we introduced the following notation 



|k + q| 



dkk / dpp 
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n{E^k))-n{EAp)) 



1 + yU/i 



k ■ (k + q) 



EM-EM 



kp 

5(p 2 ~\k + q\ 2 ). 



Integrating over 0k, we obtain 



4vr 2 



dk \ dp 



n(E ll (k))-n(E^p)) s q 2 -(k-p) 
' E,{k)-EAV) \(k + p) 2 -q 2 



2 s MM 



(17) 



(18) 



(19) 



It immediately follows from symmetry that F\-\ = The further caculations are 

different for the cases pL — fM and p = — //. In the next subsections we consider these cases 
seperately. 

C. The case p = —p 



The expression for F 1 _ 1 (g) can be conveniently rewritten in the form 

= 7^ f~ dk[n(E +1 (k))G + (k,q) + n(E-i(k))G-(k,q)], 
Ztt Jo 



where the function G±(k, q) is defined as 



G ± (k,q) 



dp 



q 2 — (k — pY 



(20) 



(21) 



lo (k + p)(k - p ± q ) \ (k + p) 2 - q 2 

It is presumed that the integrand is zero when the expression under the square root is 
negative. 

Consider now the case of small q <C kp- Then the actual limits of integration in Eq. fl2T|) 
are (k — q, k + q), i.e. p is close to k. This allows to evaluate Eq. fl2~l|) analytically 
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G±{k > q) ~^L dp k- P ± qo 



7T 



'4k 2 



qo - V<7o - <? 2 %o - g) 



(22) 



where 9{x) is a step-function. 

We see that in the small-g limit, the g and k dependencies in the G(k, q) get decoupled 
from each other so that the g-dependence of the response function F\_\{q) is the same as 
the g-dependence of Eq. (p2|) : 

n(E +1 (k))-n(E^(k)) 



Qo ~ Q 2 0(q ~ 0) 



dk- 



k 2 



(23) 



First, we realize that the response function exhibits a square root singularity at q = q . This 
singularity is depicted schematically in Fig. 2a. The second important observation is that, 
since go is determined by the parameters of the electron spectrum and does not depend on 
the Fermi level position (see Eq. fllOD), the temperature dependence of F 1 _ 1 (g) is very weak. 
Indeed, the straightforward calculation of the integral over k yields 

.2 / rp N 2l 



Fi,-i(q,T) 



Sk 2 F 



qo 



ql - q 2 6(q - q) 



7T 



1 + 



(24) 



Such temperature dependence should be contrasted to that of the conventional 2kp- 



anomaly. At finite T, the anomalous term, J q 2 — 4k 2 ?, in Eq. (Q) is replaced by [3,4] 



1 r y/g 2 -4k 2 F + 8m*Tx/h 2 

2 J cosh 2 x 



(25) 



which leads to the smearing of the anomaly within the region \q — 2kp\ ~ kpT jEp. 

Let us consider the case of large q ~ (kp^ + kp^ = 2kp. In this case the limits of the 
integration in Eq. (pl|) are (q — k, oo). By noting that the main contribution to the anomaly 
comes from the vicinity of q ~ (k + p), the integral can be simplified as 

1 r°° dp 7r 6(q - 2k =F g ) 



G ± (k,q) 



(26) 



/2g J q -k (k-p±q )^k+p-q ^2q ^q - 2k =f g ' 
Substituting Eq. ( p6|) into Eq. (p0|) leads to the following results for the response function 



v 



V2' 



for q < 2kp 



(27) 



«.-•<») =-^(1-^-1). for < 28 > 

Eqs. (p7|) and ( |28| ) apply in the domain g — 2/cp <C We conclude that the anomaly 
occurs at q — kp + kfp = 2k f. This should be contrasted to the anomalies coming from 
the terms with // = // in Eq. ( pT5|) (intrasubband transitions). These anomalies are studied 
in the next subsection. 



D. The case fj, = fj, 



In this case the anomalies occur at q = in Fi^(q) and at g = 2/^ in To 
extract the non-analytic behavior, we note that the main contribution to the integral Eq. 



(|19|) comes from k,p ~ kp and q ~ fc + p, so that we can simplify the expression for F 1: i(q) 
and F_i j _x(?) in the following way 

y/2v 



Ky/qiq + go) 



dkdp 



n(E+i(k))y/k + p - q 
k — p 



(29) 



V2i 



^Vq(q - go, 
Then by straightforward calculation we get 



dkdp 



n(E_i(k))y/k + p- q 
k — p 



F hl {q)=F-i,-M = ^ for q < 2A£' 2) , 



(30) 



(31) 
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-2A;i 1) \ 3/2 l 



2,kp ^ ' 



for q > 2k 



(i) 

F > 



(32) 



v 7 ^ 



g-24 2) \ 3/2 



2fc 



(2) 



for g > 2fc£ ) . 



(33) 



We see that instead of conventional \/q — 2kp anomaly, we obtained much weaker anomalies 

/ (12)\ 3 ^ 2 (12) 

oc I q — 2kp I . The origin of such a weakening is that for virtual transitions, kp — 2q — > 

(1 2) 

—k F ' , which give rise to the conventional Kohn anomaly, the corresponding overlap integral, 



Eq. fllq), is identically zero in the presence of the SO coupling. 



III. OSCILLATIONS OF THE SCREENING POTENTIAL 



Consider a point charge e at distance d from the 2D plane. The potential created by this 
charge within the plane is given by [3] 

As we demonstrated in the previous section, the dielectric function e(q) exhibits anomalous 
behavior at small q = qo and also at q — k$ + kp\ Besides, there are weak anomalies at 
q = 2kp\ q = 2k jP caused by the intrasubband virtual transitions. 

The general approach for extracting the oscillating behavior of V(p) is as follows. Suppose 
that e(q) has an anomally at some q = q c . Then the integrand in Eq. (|34| ) should be expanded 
with respect to bT defined as 

8F(K) = F(q c + K)-F(q c ), (35) 

and the Bessel function should be replaced by its large-p asymptotics. This leads to the 
following general formula for the oscillating part of V(p) 

ev(q c ) f~'< 



V(p) 



e V KQcP 
where the integrals ii, 7 2 are defined as 



cos(g c p - ^-)h - sin(g c p - ~)I 2 



exp(-q c d), (36) 



Ji = J dnSF^K) cos(ftp), (37) 

h = J dK5J 7 (K) sin(/tp). (38) 

By using the results obtained in the previous sections, we can now easily carry out the 
calculations (Eq. (j36l)- Eq. (j38[)). 

(i) For anomaly at q c = qo, we have 



8 k 

and e(g c ) ~ e . This leads us to the spin-orbit oscillations Eq. (|Tl 



8F{k) « (39) 



(ii) For anomaly at q c = kp + kp = 2kp, the form of bT was shown to be 

= (40) 

where we have considered the contribution of the F-i^q) which gives the same contribution 
as It is worth noting that this form coincides with the standard form of the Kohn 

anomaly in 2D. Thus for the oscillations with period n/kp we get the same form of the 
screening potential as Eq. (p. 

fl 2) _____ 

(iii) At q c = 2kp , it follows from Eqs. (|32[) and (|33|) that 

3/2 

Substituting this form into Eqs. (|36| - p8|) results in 



^w=Ki?W) 3/2 ' (41) 



3V2e 2 sm(2k ( p' 2) p) 



V °- M = (42) 

We see that oscillations with "usual" periods ir/kp' 2 ^ fall off as 1/p 3 , i.e. faster than the 
conventional oscillations Eq. (HI). 



IV. CONCLUSION 

The main result of the present paper is the prediction of the novel oscillations of the 
screening potential originating from the spin-orbit term in the Hamiltonian of a free electron. 
In our calculations we assumed that the bare energy spectrum is parabolic. We also assumed 
that the origin of the spin-orbit term is the asymmetry of the confinement potential [2] and 
neglected the term coming from the absence of the inversion symmetry in the bulk [5]. 
Note, however, that if the Dresselhaus mechanism [5] dominates the SO coupling (in this 
case Hso oc io x k x — <J y k y )), the results of the present paper remain unchanged. 

There are two cases when the spin-orbit-induced oscillations, studied in the present paper, 
dominate over the conventional 2/cp-oscillations [1]. First, with increasing temperature, 2kp- 
oscillations get damped, while go _oscma tions remain almost temperature-independent. This 
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situation is illustrated in Fig. 3. It is seen that for parameters of Fig. 3 the Friedel 
oscillations almost disappear at T ~ Q.lEp. Second, with increasing the spacing d from 
the ionized impurity to the 2D electron plane the amplitude of the conventional oscillations 
falls off as exp(—2kpd) [1], whereas for go-oscillations this dependence is much weaker oc 
exp(-qod). Fig. 2 illustrates how the go-oscillations take over as d is increased. Note, that 
large values of the spacing d is a necessary requirement for high mobility of 2D electron 
systems [6]. On the other hand, we considered the ideal electron gas, in the sense that we 
assumed the condition q l 3> 1, I being the mean free path, to be met. Thus, the large 
values d is the domain of validity of our results. 

Another application of the results obtained in the present paper is the effect of the SO 
coupling on indirect exchange interaction of localized magnetic moments (RKKY interac- 
tion). It is straightforward to generalize the expression for RKKY in two dimensions [7] to 
include the SO coupling 



Hrkky(pi2) — J 



$i(pia)Si • S 2 + ($ 2 (pi 2 ) - $i(pi 2 ))(Si • n)(S 2 • n) 



(43) 



where pi 2 is the distance between the two localized spins Si and S 2 , J is the contact inter- 
action constant, and the functions $i(p), $ 2 (p) are defined as 

*i(p) = £/^^m^)> (44) 
mm' 

* a (p) = E/705^ P ^'(ff)- (45) 
mm' 

In Eq. ( fl~5|) the response function F >(q) is defined by Eq. flilf) , whereas II '(g) is the 
conventional polarizability 

n , : ^ "(^W)-»(£y(|k + q|)) IAR , 
= E g „ W - £|/( | k + q |) ■ < 46 > 

In the absence of the SO coupling we have $i(p) = $2(p), and the interaction falls off 
at large distances as sin(2/cpp)/ p 2 . As it was demonstrated above, with SO coupling 
the function $i(p) contains high-temperature oscillations, cos(gop)/p 2 resulting from the 

11 



term Fi t -±(q) = F-i^q) in the r.h.s of Eq. fl44|) . The function $2(p) also exhibits high- 



temperature oscillations due to the term IIi^i^) = n_ 11 (g) in the r.h.s of Eq. P5j). The 
small-g anomaly in this term can be extracted in the similar way as it was done in Sec. II. 
We list only the final result 



IJ,_,(,> = 7$p4 (47) 



It is seen that the anomaly at q = q in IIii is much stronger than in F\ \. As a result, 
at large distances, $2(p) is much bigger (by a factor ~ kpp/q ) than $i(p). Then the p- 
dependence of the indirect exchange interaction takes the form sin(gop)/p- As it follows 
from Eq. (^) the in-plane components of Si and S2 interact at large p much stronger than 
the normal components. 

As a final remark, we address the above assumption that the SO splitting of the Fermi 
surface does not depend on hp- Strictly speaking, since Icf is determined by the electron 
concentration n, and, in turn, the change of n results in the change of the confinement 
potential, there should be a certain dependence of a on n. Recent magnetotransport ex- 
periments on different 2D electronic systems yield either a weakly changing with n (Refs. 
[8], [9]) or, a, roughly, inversely proportional to n (Ref. [10]). Meanwhile, it is important 
to note that a depends only on the total concentration of electrons, which does not depend 
on temperature. Thus, our main conclusion that spin-orbit-induced anomaly in e(q) is not 
smeared out by temperature is insensitive to the actual dependence a{n). This dependence 
can govern only the position of the anomaly. 
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FIGURES 

FIG. 1. The energy spectrum of a 2D system with spin-orbit coupling; upper and lower 
branches correspond to chiralities fi = — 1 and /j, = 1 respectively; the difference between two 
Fermi momenta kpP and kft is go determined by Eq. (10). 

FIG. 2. The numerical results for large-p behavior of dimensionless screened potential 
V{p){-^^)~ 1 at T = are plotted for different distances, d, of a point charge from the 2D 
plane at qo/2kp = 1/4. (a) d = as- The solid curve represents the result calculated within the 
Thomas-Fermi approximation: e(q) = eo(l + 2/qaB)- Dotted curve is the RPA result in the absence 
of the SO coupling featuring the Friedel oscillations. Dashed-dotted curve is plotted with F(q) in 
the form of Eq. (23) and features the go-oscillations. The inset shows schematically the small-g 
anomaly in the response function F(q). (b) the same as (a) for d = 1.5a#. (c) the oscillating parts 
of the screened potential for d = as and d = 1.5a#. Dotted and dashed-double dotted curves 
correspond to the Friedel oscillations. Thick solid and thick dashed-dotted curves correspond to 
the go- os ciUations. 

FIG. 3. Asymptotics of the screened potential V(p) in the units of 2e 2 /e^as are plotted at 
different temperatures at qo/2kp = 1/4. Dotted and solid curves show the 2/c^-oscillations at 
T = and T = O.lEp respectively. Thick curve represents the go-oscillations, Eq. (11), their 
change with temperature from T = to T = O.lEp being smaller than the thickness of the curve. 
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